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ABSTRACT 

The present paper extends the method of computing the natural modes of 
flexural vibration of a nonrotating propeller blade which was described in RP764 
to blades that rotate as in actual flight. The resulting integral equations were 
solved for two aluminum-alloy propeller blades of typical design vibrating with 
the fundamental mode and with the second harmonic mode and rotating at speeds 
covering the range of service speeds. The effect of rotation on the stress distri- 
butions and on the natural frequencies was obtained for the two extreme end 
conditions of rigid clamping at the hub and no clamping at the hub. Rotation 
was found to have a small effect (up to 13 percent) on the maximum stress per 
unit tip deflection and was found to shift the maximum toward the hub. For 
the determination of critical speeds, the natural frequencies were shown to follow 
conveniently from Lord Rayleigh's energy method by making use of the solutions 
for no rotation outlined in RP764. The effect of initial twist on the flexural 
frequencies of the blade was found to be small compared with the effect of changes 
in clamping at the hub. Changes in clamping may change the frequency up to 
13 percent. A correction to RP764 is noted in this connection. Both the stress 
distribution and the natural frequencies of a given propeller blade may be calcu- 
lated from those obtained for a model blade by means of two model rules developed 
by Theodorsen. 
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I. PURPOSE OF INVESTIGATION 

The determination of stresses in a nonrotating propeller blade 
vibrating with a flexural mode has been discussed in Bureau Research 
Paper RP764 [I]. 1 In that paper it was shown that the bending 
stresses set up in the blade could be conveniently determined from 
strain measurements on the surface of the blade, after the blade had 
been set into resonance with the desired flexural mode. The stresses 
so determined were found to agree closely with those calculated from 
the theory of flexural vibrations for beams of variable section. Fatigue 
failures in the center portion of eight aluminum-alloy blades were 
produced by vibrating the blades at the fundamental mode with a 
sufficiently high amplitude to produce a crack after 33,800 to 2,390,000 
cycles of vibration. The crack occurred in every case at a maximum 
stress amplitude greater than 14,000 pounds per square inch, and at 
a point where the stress differed only a few percent from the measured 
maximum stress. Attempts to produce fatigue failures near the tip 
in blades vibrating with the second harmonic mode were not success- 
ful, mainly because a sufficiently high stress amplitude could not be 
maintained in the blade. Within the last few months one such fatigue 
failure was obtained in an aluminum-alloy blade while it was being 
used to calibrate certain dynamic-strain gages. The fatigue crack 
extended through about 85 percent of the cross-sectional area of the 
blade and developed at a distance of 13 inches from the tip at a point 
within an inch of the location at which the maximum stress amplitude 
had been measured. 

The results of tests on nonrotating blades can be applied to propel- 
lers in flight only if the effect of centrifugal force on the natural modes 
of vibration is clearly understood. A number of investigations have 
already been made to determined the effect of an angular speed, o>, 
on the natural frequency, p, of the blade, and a formula of the type 

p*=p Q 2 +au 2 (md 2 /sec 2 ) (1) 

was proposed as early as 1918 by Berry [2]. Values of the constant 
a ranging from 1.04 to 1.9 have been obtained by various investigators 
[3, 4, 5, 6, 7, 8, 9]. The large variation in the values of a may be 
ascribed, principally, to differences in the design of the propeller 
blade. The constant a was determined theoretically as 1.19 for the 
uniform beam vibrating with its fundamental mode, by Berry [2], 
and as 1.4 for a wedge, by Webb and Swain [3]. Hansen and Mesmer 
[7] concluded, on the basis of tests, that a is about 1.45 for blades of 
usual design. Theodorsen obtained a considerably higher value 
(<*=1.9), also on the basis of tests (revision of values given in [9].) 

The effect of centrifugal force on the stress distribution in a rotat- 
ing blade vibrating with a flexural mode has not been studied as 
extensively as the effect on the natural frequencies. An experimen- 
tal determination of the stress distribution has not been possible 
until very recently, since no dynamic strain gages were available 
which could be attached to a rotating blade and which could be 
expected to give a fairly accurate value of the strain amplitude dur- 
ing lateral vibrations. In addition, the problem of exciting a single 

1 Figures in brackets indicate the literature references at the end of this paper. 
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flexural mode in a rotating blade is not a simple one; up to the pres- 
ent time it has been solved only for a thin model blade rotating at a 
relatively low speed [91 . 

Until the experimental difficulties just mentioned have been over- 
come, an answer to the problem can only be expected from a theoret- 
ical attack. The close agreement between measured and calculated 
stress distributions and natural frequencies for nonrotating blades 
suggested that an extension of the solution, to rotating blades would 
give equally reliable results on the effect of rotation on both stress 
distribution and natural frequencies. 

The present paper deals with the extension of the theoretical method 
originally developed by Hansen and Mesmer [7, 1] to rotating blades 
and gives the resulting stress distributions and natural frequencies for 
two aluminum-alloy propeller blades of different design. 



II. CALCULATION OF STRESSES 

1. EQUATION OF MOTION OF IDEALIZED STRAIGHT PROPELLER 

BLADE 

It will be assumed, as in the case of the nonrotating blade [1], that 
the flexural vibrations of the propeller blade may be adequately rep- 
resented by the motion of a straight beam whose cross-sectional areas 
and principal moments of inertia are equal to those of the propeller 
blade, and whose principal axes of inertia remain parallel throughout 
the length of the beam. The ordinary beam theory may then be 
applied in setting up the equation of motion of the blade for the case 
of resonant vibration in the fundamental mode and the second har- 
monic mode. 

The adequacy of this approximation as applied to nonrotating pro- 
peller blades was confirmed by the close agreement that was found [1] 
between the stress distribution and the frequency of the actual blade 
with its initial twist and initial curvature and the corresponding stress 
distributions and frequencies for the idealized straight blade. Its 
adequacy for determining the fundamental frequency and second har- 
monic frequency of rotating blades is confirmed by the close agreement 
between the experimental frequencies determined by Theodorsen (re- 
vision of values given in [9]) and the frequencies calculated assuming 
an idealized straight propeller blade (see fig. 17, page 660). 

No data are available to verify the assumption that the stress dis- 
tribution for the fundamental mode and the second harmonic mode 
of the rotating blade is not affected by the initial twist. Strain meas- 
urements to check this assumption are desirable. 

A small section of the idealized straight blade of length Ax (fig. 1) 
will be in equilibrium in the ^-direction under the action of two 
forces; (1) the inertia force pA'zAx, where p is the density and A is 
the cross-sectional area of the blade; and (2) the restoring force in 
the 2-direction tending to bring the blade back to its mean position 
due to the tension, T, and shear, Q. This condition may be expressed 
mathematically by 

g- (T sin 6+Q cos 6) Ax=pA'zAx. (2) 
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For small deflections, only terms of the first order in need be re- 
tained, so that sin d^'bzl&x, cos — 1 and therefore eq 2 may be 
written as 

The ordinary beam theory gives with the same degree of approxi- 
mation 



d d 2 z 



(4) 



The tension, T, at x due to centrifugal force is for small deflections 
(cos 0—1) given by the sum of the tensions pAxAxco 2 due to the par- 




Figure 1. — Forces on portion of rotating propeller blade. 

tides lying between section x and the tip of the blade. This leads to 
the integral 

T=pa> 2 ( L Axdx. (5) 



v*jr 



Substituting eq 4 and 5 in eq 3, and dividing through by Az, gives the 
following differential equation for the motion of a rotating beam 
vibrating flexurally 



dx 



[-i( £/ S)+^ir^]"^ =s °- {6) 



During a natural vibration of a blade, each of its sections will 
vibrate sinusoidally with the natural frequency, p, of the mode con- 
sidered. If X is the amplitude of motion of the section, its motion may 
be described by 

z=X sin pt. (7) 
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The substitution of eq 7 in eq 6 leads to the following differential 
equation for the natural modes of flexural vibration of a rotating blade: 



d 

dx\ 



[-^ i w)+^XH +p ^ x=0 - (8) 



It is convenient in considering solutions of this equation to introduce 
the following dimensionless variables and the constant k: 

„ x X r A • I j pL 2 p' 

*-Z' ' = Z' 



H? i= h k=e if' «=//*#• < 9 > 



The ratios are proportional as follows: £, to the distance from the 
axis of rotation of the section considered; rj, to the amplitude at the 
section ; /, to the cross-sectional area ; i, to the moment of inertia ; Jc, 
to the square of the natural frequency; and u, to the centrifugal ten- 
sion. In terms of these variables, eq 8, after a simple transformation, 
becomes 

§^=</«+£|(W)) do) 

where the prime denotes differentiation with respect to £. 
The amplitude of the bending moment at any section, 

M=EJ^^ELHn', (11) 

may be derived from this equation by integrating twice with respect 
to £. The two constants of integration are fixed by the boundary 
conditions of zero shear and zero bending moment at the tip. These 
are, 

e-if-^wo-OpV-o. (12) 

Integrating eq 10, and substituting the boundary conditions, in eq 12, 
gives 

v "4J* I J» ljM& *~?X**'*} (13) 

The first integral on the right denotes the bending moment due to 
inertia forces, while the second integral is the restoring moment due 
to the centrifugal forces. This latter restoring moment is found to 
increase with the square of the speed of rotation of the blade. It 
adds materially to the stiffness of the blade at speeds of rotation 
approaching operating speeds. 

The extreme fiber bending stress, a x , is according to the ordinary 
beam theory, 

v x =-Y=jEr\\ (14) 

where c is the distance of the fiber from the neutral axis of the section, 
and 17" is as given by eq 13. 



- 
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An integral equation for the amplitude curve r?(£) can be obtained 
by integrating eq 13 twice again and using the proper boundary 
conditions at the hub. These boundary conditions will depend upon 
the type of exciting forces and restraints. Because of the twist in 
propeller blades, the natural modes are neither purely symmetric 
nor purely antisymmetric. However, the motion for the lower-fre- 
quency modes is mainly parallel to the shaft, and the types of motion 
may be classified according to the symmetry or antisymmetry of this 
component. 

Modes excited by forces at the center of the hub in the plane of 
rotation and modes excited by moments about axes in this plane will 
have components normal to this plane which are antisymmetric, 
while modes excited by forces at the hub parallel to the shaft and 
modes excited by moments about the shaft will have components 
normal to the plane of rotation, w r hich are symmetric. Since the 
principal source of vibration in propellers is believed to be torque 
variation in the shaft, only the second type of motion is considered 
important. On this account, only symmetric modes in the idealized 
straight propeller are considered in this paper. 

Integrating eq 13 twice for a symmetrical mode of vibration, that 
is, a mode of vibration satisfying the end condition, 

£=0, 9 '=0, (15) 

gives the equation 

™+ k E$[f!P*^£**)?* (16) 

The constant of integration r) , which corresponds to the deflection 
at the hub, will depend on the restraint at the hub. The hub restraint 
at the root of each blade may be approximated by setting it equal to 
that due to a concentrated mass, M, rigidly attached to the hub com- 
bined with an elastic restraint against axial motion with a spring 
constant, S. The harmonic motion with an amplitude X =Lrj of 
the mass, M, will apply a deflecting force with an amplitude MXpp 2 
to the hub, while the elastic restraint applies a restoring force with 
an amplitude SX . Balancing the resultant of these forces by the 
shear in the hub section of the blade gives 

JW-flQt-g(jSlg)]^ (17) 

It is convenient to express M and S in terms of the ratios, 

M S n^ 

m== W s= ^f' (18) 

so that eq 17 becomes with eq 9, 

^-•>-G¥L (19) 

It will be noticed from eq 13 that for a symmetrical mode, ??'(0)=0, 

H3L~J> 
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so that eq 19 becomes 

7j (m-s) = - ( fyd£ (21) 

This equation evidently expresses the equilibrium between the con- 
centrated transverse force applied at the hub and the distributed 
transverse inertia forces due to the vibration of the blade integrated 
from the hub to the tip. Substitution from eq 16 for the right-hand 
side of eq 21 leads to an equation which may be solved for r} . Sub- 
stituting the resulting value of i] in eq 16 leads to the following 
integral equation for symmetrical flexural modes of vibration: 

Numerical solutions of this integral equation were obtained for two 
special cases only; (1), the case of rigid clamping (fixed-free blade), 
which corresponds to infinite values of m— s } so that eq 22 becomes 



i=k 



x^[n>^ (2a) 

and (2), the case of no clamping (free-free blade), which corresponds 
to ra— s=0,*so that eq 22 becomes 

It is interesting to note that the boundary condition m— s=Q for a 
free-free blade may be satisfied theoretically in the presence of elastic 
restraint at the hub, provided the concentrated mass at the hub is 
such as to balance the spring force at the natural frequency in question. 

The solutions, eq 23 and 24, will bound all those natural modes for 
which m— s>:0, that is, they apply to the symmetrical flexural "modes 
in a two-blade propeller whose hub acts like a concentrated mass, 
Mq=pI? (?n—s). In the remaining cases in which the propeller hub 
acts like a spring with an equivalent spring constant, S =pL 3 p 2 (s—m) f 
the frequency will, in general, lie outside of the values determined by 
eq 23 and 24. If, for example, in eq 22 m-s is given a large but 
finite negative value, the solution will have a frequency below that 
determined by eq 23; while if m—s is given a small negative value, 
the frequency will be above that determined by eq 24. 

Experiments indicate that the natural flexural modes of non- 
rotating propellers [1} are bounded by the particular solutions in eq 23 
and 24 with co=0. Propeller hubs of practical design seem to act 

100926—38 8 
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like a concentrated mass rather than like a concentrated spring. It 
was decided, accordingly, to confine the theoretical work on the 
rotating blade to solutions of eq 23 and 24 for certain propeller blades 
of practical design. The resulting modes were found to be sufficiently 
similar for the two boundary conditions of rigid clamping and no 
clamping at the hub, as expressed by eq 23 and 24, to make it prob- 
able that they are adequate expressions for the flexural modes of the 
blades in a multiple-blade propeller as well as in a two-blade propeller. 

2. MODE OF MOTION AND STRESS DISTRIBUTION FOR 
FUNDAMENTAL MODE 

Equations 23 and 24 may be solved for the fundamental mode, and 
the stress distribution for this mode, by using a procedure that is a 
direct extension of the corresponding procedure for nonrotating 
propeller blades, as described in detail in Bureau Research Paper 
RP764 [1]. The deflection curve r? for the fundamental mode for a 
given value of the ratio o)/p is obtained by assuming a deflection 
curve rj which approximates roughly the shape of the expected 
deflection curve, substituting the assumed value of t\ on the right- 
side of eq 23 or 24 and carrying out the integration by the numerical 
procedure described in detail in [1]. This will, in general, lead to a 
value of rj/k which approximates the shape of rj for the fundamental 
mode more closely than the initially assumed shape of rj. The 
integration is repeated after inserting on the right-hand side of the 
equation a deflection curve proportional to the deflection curve rj/k 
just calculated, until the integral r\\k has the same shape as the r? in 
the integrand on the right. The constant of proportionality, k, will 
then give the angular frequency, p, of the blade, and knowing p, 
one may calculate the speed of rotation, co, from the initially assumed 
value of u/p. The calculation was found to converge to the correct 
solution after two or three cycles of integration when the initial rj 
was taken as the deflection curve of the propeller blade without 
rotation. With the calculation of the deflection curve for the funda- 
mental mode of vibration the stress distribution may be determined 
by substituting the integral, eq 13, which is evaluated automatically 
as one of the steps in the integration of eq 23 or 24, in the stress 
formula, eq 14. 

It should be noted from eq 23 and 24 that u/p must be below a certain 
value in order to lead to a positive value of k and hence a real natural 
frequency, p. This limiting value of ujp corresponds to the asymptote 
to the curve ofp=f(oo) for infinitely high speeds of rotation, co. In pro- 
peller blades of usual design the asymptotic portion of the curve lies 
far beyond the range of operating speeds; the limiting value of co/p is, 
therefore, not of practical importance in this connection. 

The above procedure was used to determine fundamental modes of 
vibration and stress distributions corresponding to these modes for 
aluminum-alloy blades of two typical designs. One design, which is 
referred to as type 32, had already been used to study flexural vibra- 
tions of a nonrotating blade [1]. The other design will be referred to 
as type 4371. A model blade of this design was used by Theodorsen 
[91 to determine the natural frequencies of flexural vibration with 
rotation. The radius, L (distance from center of hub to tip), for a 
propeller blade of type 32 was equal to 47.6 inches, while it was equal to 
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66 inches for a blade of type 4371. The cross-sectional area, A, the 
moment of inertia, /, the centrifugal-force ratio, u, and the extreme 
fiber distance, c, for the convex side for any section were evaluated 
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Figure 2. — Distribution of (A/L 2 ) 1 ^ along propeller blades of type 82 and type 
4871, where A— cross-sectional area at distance xfrom center of hub and L= length 
of blade. 

from the design data for the two types of blades and were used to plot 
the curves of JM, iH { u y *, and c/L shown in figures 2 to 5. Figure 6 shows 
the ratio u/f, which is proportional to the stress due to centrifugal 
force. 
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Figure 3. — Distribution of (I/L 4 ) 1 '* along propeller blades of type 82 and type 4871, 
where I ' = minimum moment of inertia at distance x from center of hub and L = 
length of blade. 

Figure 7 gives the shape of the fundamental mode for a blade of 
type 32 calculated from eq 23 for speeds between and 2,250 rpm, 
assuming clamping at the center of the hub. The effect of centrifugal 
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force is clearly seen in reducing the curvature in the tip portion of the 
blade and in shifting the point of maximum curvature toward the 
hub. 
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Figure 4 
where u 
center of hub, A = cross-sectional area, and L— length of blade 



Figure 8 shows the corresponding result for a blade of type 4371 for 
speeds of rotation of 0, 860, and 2,020 rpm. 
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Figure 5. — Distribution of c/L along propeller blades of type 82 and type J+871, 
where c= extreme fiber distance for convex side along center line of blade at a distance 
x from center of hub and L — length of blade. 

Figure 9 shows the shape of the fundamental mode for a blade of 
type 32, as calculated from eq 24, for the condition of zero fixity at the 
hub for speeds of rotation of and of 1,620 rpm. Rotation has a very 
small effect on the shape of the mode in this case. 
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The effect of rotation on the shape of the stress-distribution curve 
for the fundamental mode is brought out in figures 10 to 12. 

Figure 10 gives the stress distribution along the center line on the 
convex side for a blade of type 32 clamped at the center of the hub for 
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Figure 6. — Distribution of centrifugal stress ratio u/f along propeller blades of type 
32 and type 4871; for definition of u and f> see figures If. and 2; L— length of blade. 

speeds of rotation between and 2,250 rpm. The point of maximum 
stress amplitude is shifted towards the hub by the stiffening action of 
the centrifugal force. At the same time, the maximum stress per inch 
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Figure 7. — Deflection curves for fundamental mode of blade type 82 clamped at 
center of hub, rotating at a speed N, and vibrating with a frequency f. 

Curve a, N*=$ rpm, /=30.6 c/s; curve 6, iV- 1,060 rpm,/=39.2 e/s; curve c, N- 1,550 rpm,/=47.1 c/s; curve d, 
iV=l,960 rpm, /=54.6 c/s; and curve e, A r =2,250 rpm, /=60.2 c/s. 

tip deflection is increased somewhat to compensate for the lessening 
in curvature and hence in stress in the outboard portion of the blade. 
Figure 11 shows the effect of rotation on the stress distribution in 
the same propeller blade, assuming the hub to bejiree. 
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Figure 8. — Deflection curves for fundamental mode of blade type 1+371 clamped at 
center of hub, rotating at a speed N, and vibrating with a frequency /. 

Curve a, N=0 rpm, /=26.0 c/s; curve 6, iV=860 rpm, /=32.6 c/s; and curve c, iV=2,020 rpm, /=50.5 c/s. 

Figure 12 shows the effect of rotation on the stress distribution of a 
blade of type 4371 clamped at the center of the hub. The shape of 
this blade is such as to lead to two points of maximum stress amplitude. 
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Figure 9. — Deflection curves per inch tip deflection for fundamental mode of blade 
of type 32 rotating at a speed N and vibrating with a frequency /, but without clamp- 
ing at center of hub. 

Curve a, N=0 rpm, /=34.1 c/s; and curve 6, iV=l,620 rpm, /=55.2 c/s. 

Deflection curve per inch tip deflection for nonrotating blade with rigid clamping at 

center of hub. 

Curve c, 2V=0 rpm,/=30.7 c/s. 

The first maximum lies at a section 55 percent of the blade length 
from the hub and predominates at low speeds of rotation, while the 
second maximum, which is 35 percent of the blade length from the 
hub, predominates at high speeds of rotation. In spite of a consider- 
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Figure 10. — Extreme fiber stress per inch of tip deflection along the center line on the 
convex side for the fundamental mode of blade type 32 clamped at center of hub, 
rotating at a speed N f and vibrating with a frequency f. 

Curve a, 2V=0 rpm, /= 30.6 c/s; curve b, N= 1,050 rpm,/= 39.2 c/s; curve c, N= 1,550 rpm,/=47.1 c/s; curve d, 
iV=l,900 rpm, /=54.6 c/s; and curve e, iV=2,250 rpm, /=60.2 c/s 
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Figure 11. — Stress distribution for fundamental mode of blade type 32 rotating at a 
speed N and vibrating with a frequency /, but without clamping at center of hub. 

Curve a, N=0 rmp,/=34.1 c/s; and curve b, iV=l,620 rmp, /=55.2 c/s. 

Corresponding curve for nonrotating blade with rigid clamping. 

Curve c , iV=0 rpm, /=30.7 c/s. 
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able change in the shape of the stress-distribution curve, the magni- 
tude of the maximum stress per unit tip deflection is changed less 
than 13 percent as a result of rotation. 

Figures 10, 11, and 12 all lead to the same conclusion, namely, that 
the effect of speeds of rotation in the operating range (1,000 to 2,000 
rpm) on the maximum stress per unit tip deflection for the funda- 
mental flexural mode may be neglected when designing propellers, 
especially in view of the uncertainty in the value of the tip deflection. 
It is allowable for approximate calculations to take the resultant 
stress in the blade as the sum of the static centrifugal-force stress 
(fig. 6), the static bending stress due to air forces, and a sinusoidally 
varying stress whose amplitude may be calculated, provided the tip 
amplitude of the blade is given, by using the same ratio of stress to 
tip amplitude as that determined either by direct measurement or 
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Figure 12. — Extreme fiber stress per inch tip deflection along the center line on the 
convex side for the fundamental mode of blade type 4371 clamped at center of hub, 
rotating at a speed N, and vibrating with a frequency f. 

Curve a, iV=0 rpm, /= 26.0 c/s; curve 6, iV=860rpm,/=32.6 c/s; and curve c, iV=2,020 rpm, /=50.5 c/s. 

by calculation, upon the same blade vibrating with the fundamental 
flexural mode without rotation. 

The preceding statement assumes that the steady air forces have 
no effect on the vibration of the propeller blade. This assumption 
is supported by the theoretical work of Reissner [17], who found that 
the modes of rotating propeller blades are practically unaffected by 
small initial curvatures such as may be produced by air forces. 

If the location of the point of maximum stress is desired the possi- 
bility of a considerable shift toward the hub must be considered in 
those cases (e. g., fig. 12) in which the stress for no rotation decreases 
slowly on the hub side of the point of maximum stress. 

3. MODE OF MOTION AND STRESS DISTRIBUTION FOR SECOND 

HARMONIC MODE 

(a) OUTLINE OF PROCEDURE 

The procedure outlined above for calculating the fundamental 
mode cannot be directly applied to the calculation of one of the har- 
monic modes. The procedure is a method of iteration in which the 
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solution generally converges to the fundamental mode. The diffi- 
culty may be circumvented, at least for the second harmonic mode 
[1], by forcing the deflection curve through a nodal point and vary- 
ing the position of the nodal point until eq 23 or 24 is satisfied for all 
sections along the vibrating blade. Obviously, this is a very cumber- 
some procedure and one that can hardly be extended to a calculation 
of harmonics with more than one nodal point along the blade. 

A more convenient and more generally applicable procedure may 
be derived by eliminating the modes of lower frequency than the de- 
sired mode from the assumed mode. The procedure is easily under- 
stood if it is realized that the solutions of eq 22 are orthogonal. This 
may be proved as follows. 

Let <p p (£) t cp Q (0 be two solutions -q of eq 22. They will then satisfy 
eq 10 independently, and we may write eq 10, first with <p p in place of 
?7, and then with <p q in place of tj. Multiplying the first one of these 
equations by <p q and the second one by <p P) subtracting one equation 
from the other, and integrating over the length of the blade gives the 
following equation: 






(25) 



Integrating the terms on the right side of this equation by parts leads 
to 

-CH^toN*'-*!*/) • (26) 

The terms on the right of eq 26 disappear for the upper limit £=1, 
since the shear force (proportional to di<p"/d£), the bending moment 
(proportional to i<p"), and the centrifugal force (proportional to u) 
must be zero at the tip f=l. At the lower limit, |=0, all terms except 
the first pair disappear because the slope (proportional to <p') is zero 
at the hub £=0 for symmetrical modes. The first pair of terms is from 
eq 18 and 19 

= — <P Q (0)(p p (0)[m(k p —k Q )—k p s p +k Q s q ]. (27) 
From eq 18 and 19 we note that 

fcs=y-p==constant, (28) 
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so that the last two terms in eq 27 drop out. Substituting the result- 
ing value of eq 27 in eq 26 gives the relation of orthogonality 

m(p J tf> f )*„ + l/^«#=0, p^%. (29) 

The first step in the procedure of solving eq 23 and 24 for harmonic 
modes consists in obtaining an assumed deflection curve y which satis- 
fies the four boundary conditions, eq 12, 15, and 21. It is possible to 
express this curve t\ in terms of the normalized deflection curves 2 
?i({), . . • Vnii) • • • for the natural modes of the propeller blade 
that are solutions of eq 22: 

V®=a 1 < fl (0+ . . . +a n <p n (0+ .... (30) 

If the right-hand side of eq 30 is substituted in eq 22 and the integra- 
tion performed, a new curve, 

1 >= a %+ a % +a %+ ■-' (31) 

is obtained, in which the lower orders are more predominant than in 
the original curve rj, since the frequency ratio, k, increases rapidly 
with the order of the mode. The second harmonic frequency, for in- 
stance, is roughly three times the fundamental frequency for a typical 
fropeller blade; this leads to a value k 2 about 10 times as large k%. 
t is thus understood how, by reiterated integration, t) converges to 
the mode having the lowest frequency. 

A general method eliminating all modes below the one to be found 
consists in evaluating the tip amplitudes a* f « • • , a ni • • • of the 
component modes by first multiplying eq 30 by/^i(£), • • • f<p n (£), 
• • • , integrating, and then making use of the condition of ortho- 
gonality, eq 29. This leads to the values 



ar- 



a n z 



mu(0)*i(0)+ f l M(&v(&di 

Jo 

Jo 
• (32) 

mij(0)*,(0)+ f/w,tt)ntt)<# 

jo 



fj<Pn 2 WZ 



These equations enable us to free the assumed deflection curve, 
77(f), of am^ component modes whose shape has already been deter- 
mined. If we know all modes below the n-ih mode we may calculate 
the n-th mode by starting out with a deflection curve from which all 
lower modes have been eliminated, 

Vi=y—ai<Pi ■ • • — an-npn-i—anVn+an+ipn+i • • • , (33) 



2 A deflection curve, ?i (£),.. . <p n (£), will be considered as normalized in this discussion if its value at 
the propeller tip is taken as unity, that is, <p x (1)=1, . . . *>»(1)=1. 
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and substituting this deflection curve in the integrand on the right- 
hand side of eq 22, 23, or 24. Carrying out the integration leads to a 
deflection curve, 

I,=a n ^+a n+1 ^+ • • • , (34) 

in which the n-ih mode is more predominant than in 77,, since the fre- 
quency ratio, k, increases rapidly with the order of the mode. 

To get a better approximation of the n-th mode, one may normalize 
Ij by multiplying it by a factor which will make the tip deflection 
unity, and may then repeat the integration with the resulting deflection 
curve, rjjj, 

*,ft)=I,({)/I,(l) (35) 

until the resulting integral is proportional to the rj(g) in the integrand 
within the accuracy of the numerical procedure used. 

A second method of calculating the harmonic modes is considerably 
less general than the first method, but possibly more convenient, if 
only the fundamental need be eliminated. 

In the second method we free the assumed deflection curve r?({) of 
any component modes whose shape has already been determined 
after the integration of eq 31. The amplitudes a 3 , a 4 , . . . will be 
small compared to a 2l provided that 77 has approximately the shape 
of the second harmonic mode. If we let a 3 =a 4 = . . . =0, we may 
solve eq 31 for k 2 as follows: 



h = 



a<2<P2 __Vr—ai<Pi 



T Oitpi T Oxtpi (36) 

The frequency ratio, k 2 , must have the same value for all sections 
along the blade. All quantities in eq 36 are already known except 
&i and k 2 . The amplitude a x may, therefore, be evaluated by de- 
riving eq 36 for any two sections of the blade. A number of trials 
showed that the sections f— 0.45, £=0.85 were particularly suitable. 
That is, a x was evaluated from 

77,(0 .45) -a iyi (0.45) = 17/(0-85) -6^(0.85) 

JV(0.45)-^f' 45) 7,(0.85) - a ^f- 85) ' (37) 

Solving eq 37 for a x gave 

,(0.45M0.85)M-M 

n L*7(0-85) ?7(0.45)J /oo^ 

1 , 1 (0.45)[/(0.85)-^]-, 1 (0.85)[/(0.45)-^p] ^ 

Knowing a u a further approximation to tj was obtained from eq 30 by 
evaluating 

Ii~<hf (39) 

and normalizing this expression. Repeating with the normalized 
curve rj u proportional to eq 39 gave a still closer approximation t\ lu . 
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It was found that the fundamental component a x was negligible 
compared with a 2 after the third run; three runs generally sufficed to 
obtain an accurate solution of eq 23 or 24 for the second harmonic 
mode. 

A check determination of the second harmonic mode of a blade 
type 4371 rotating at 860 rpm was made to verify that both the 
procedures outlined above led to the same numerical result. 

(b) RESULTS 

The second harmonic mode was calculated by the second procedure 
outlined above for a number of rates of rotation for a blade of type 32 
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Figuee 13. — Deflection curves for second harmonic mode oj blade type 32 clamped 
at center of hub, rotating at a speed N, and vibrating with a frequency /. 

Curve a, N=0 rpm, /= 119 c/s; curve b, iV=814 rpm, f=l24 c/s; curve c, 2V= 1,210 rpm, /= 131 c/s; and curve 

d, JV=2,240 rpm, /= 153 c/s. 

clamped at the center of the hub. The results are shown in figures 
13 and 14. 

Figure 13 gives the shape of the mode for rates of rotation from 
to 2,240 rpm. The deflection is seen to change only slightly because 
of rotation. 

The effect of the rotation on the stress distribution is shown in 
figure 14. Comparison of this figure with figure 10 shows that the 
effect of rotation on the stress distribution is less than for the funda- 
mental mode. 

The effect of rotation on the natural modes may, in general, be 
expected to decrease with increased orders of vibration. In the 
higher orders practically all the energy stored in the deflected posi- 
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tion is energy of bending and the proportion of restoring force due 
to rotation becomes small compared with that due to bending. It 
may be concluded that except for the fundamental mode the stress 
distribution for a given tip deflection of a rotating blade vibrating 
with a flexural mode may be approximated closely by superimposing 
the vibratory stress in that mode, as determined for no rotation on 
the static stresses due to centrifugal force and to air loads. 
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Figure 14. — Extreme fiber stress per inch tip deflection along the center line on the 
convex side for the second harmonic mode of blade type 82 clamped at center of 
hub, rotating at a speed A T , and vibrating with a frequency f. 

Curve a, N=0 rpm;/=H9 c/s; curve b, iV=814 rpm,/=124 c/s; curve c, N= 1,210 rpm,/=13l c/s; and carve 

d, iV"=2,240 rpm, /=153 c/s. 

III. CALCULATION OF NATURAL FREQUENCIES 

An accurate knowledge of the natural frequencies of vibration of a 
propeller blade is required whenever there is a possibility of resonance 
between the propeller vibration and a periodic exciting force of suffi- 
cient amplitude from the engine. Accurate values of the frequency 
for the extreme end conditions of rigid clamping and of no clamping 
at the hub can be obtained by solving the integral eq 23 and 24, as 
was shown above. This involves too much calculation for practical 
design and fortunately is not necessary for speeds of rotation of a 
practical order. 

1. APPROXIMATE CALCULATION BY RAYLEIGH'S METHOD 

For speeds of rotation of a practical order, the shape of the deflection 
curve is not changed appreciably by the centrifugal force, and an 
accurate value of the natural frequency may be obtained, as noted 
below, by the Rayleigh energy method [10]. 

Rayleigh's method determines the natural frequency for a given 
mode X(x) of an elastic system from the equality of the potential 
energy, K, stored in the system at the extreme position and the kinetic 
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energy, T, acquired by the system when it swings through its mean 
position. The potential energy of a rotating straight bar of variable 
section vibrating transversely consists of two terms [12]: (1), the 
potential energy due to bending, 

v, =l!XW^ (40 > 

and (2), the potential energy V" due to the inward motion of the 
particles in the blade against centrifugal force. 

The kinetic energy for sinusoidal vibrations, as given by eq 7, will be 

T=^ CpAXWx. (42) 

The integrals on the right-hand side of eq 40 to 42 may be written 
more conveniently in terms of the ratios in eq 9. Doing this, and 
substituting the resulting expressions in the condition of conservation 
of energy, 

V' + V" = T, (43) 

gives the following values for the square of the natural frequency, p, 
in radians per second: 



2 P# 
p *== £ — 



E 2 \\" 2 dk+^^^\'H^di (44) 



Equation 44 may be written in the form of the well-known eq 1, 

p 2 =pi+ao) 2 (r Sid/sec) 2 } (45) 



where 



pg=£L2o ;«=ALJo — J_ 



fa*** ffa* 



provided 7j is taken as proportional to the deflection curve for no 
rotation. 

The quantity p 2 will be the square of the frequency without rotation 
and the ratio a will be constant as long as the shape t?" (£)A?(1) of the 
curvature curve, the shape r}'(£)lrj(l) of the slope curve, and the shape 
ri (£)/?7 (1) of the deflection curve remain unaffected by the rotation co. 
The detailed analysis of the previous section showed that thisisapprox- 
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imately the case for propellers of usual design rotating at service 
speeds. t . 

The closeness with which Rayleigh's method will give the correct 
solution of eq 23 and 24 for the fundamental frequency of a pro- 
peller blade, provided the mode for no rotation is substituted for rj, is 
brought out by a comparison of curves (a) and (b) in figure 15 (A) and 
15 (B) with the open circles which were calculated for a blade of type 
32 by direct solution of eq 23 and 24. 

The insensitiveness of Rayleigh's method to variations in the shape 
of the deflection curve is so great that the deflection curve, n, for a 
vibrating beam is frequently approximated by the deflection curve of 
the beam under its own weight. The deflection curve of the propeller 
blade as a cantilever beam deflected by its own weight follows from 
eq 23 by letting co=0 and by replacing the vibrational acceleration, 
p 2 X=Lp 2 rj, by the acceleration of gravity, g. Curve c in figures 




70 40 

n, rev/sec 



600 /60O 

n* (re\//6ec) 2 



Figure 15. — Variation of natural frequency, f, of fundamental mode with speed of 
rotation, n, for blade type 82. 

The points were calculated by solving integral equations 23 and 24, respectively. 

Curve a, Rayleigh's method using the mode for no rotation of a fixed-free blade, / 2 =30.9 2 -fl.95n 2 (c/s) 2 ; 
curve 6, Rayleigh's method using the mode for no rotation of a free-free blade, / 2 = 34. l 2 -r-2.46w 2 (c/s) 2 ; and 
curve c, Rayleigh's method using the deflection curve of the blade under its own weight, /" 2 =32.7 2 -}-2.04n 2 
(c/s) 2 . 

15 (A) and 15 (B) results for a blade of type 32. It is seen that the 
error in frequency as given by this procedure may be as much as 
5 percent in this particular case. 

The exactness with which the linear relation eq 45, between p 2 = (2irf) 2 
and u 2 =(27rn) 2 holds is brought out in the plot of/ 2 against n 2 in 
figure 15 (B). 

It may be concluded that eq 45 will give an excellent approximation 
to the natural frequencies of a rotating blade. for ordinary values of w, 
provided p Q is either directly observed or is calculated by the relatively 
exact method [1, 7] for no rotation. The ratio a may then be calcu- 
lated by substituting the deflection curve, ??(£), for no rotation in 
eq 47. It may be expected to vary with a change in the design of the 
blade which would change the shape of the deflection curve, n(£), or 
the weight distribution, /. 
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2. RESULTS FOR FUNDAMENTAL FREQUENCY AND SECOND HAR- 
MONIC FREQUENCY 

Rayleigh's method, using the shape of the natural mode for no 
rotation, was applied to calculate the fundamental frequency and the 
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-Variation of natural frequency, f, of second harmonic mode with speed 
of rotation, n, for blade type 82. 

The points were calculated by solving integral equations 23 and 24, respectively. Curve a, Rayleigh's 
method using the mode for no rotation of a fixed-free blade, p= U7.5H-6.8ft 1 (c/s)'; and curve 6, Rayleigh's 
method using the mode for no rotation of a free-free blade,/ 2 = 13O J +S.0n 2 (c/s) 1 . 
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Figure 17. — Variation of fundamen- 
tal and second harmonic frequencies, 
f, with speed of rotation, n, for blade 
type 4371. 

The open points were calculated by solving in- 
tegral equations 23 and 24, respectively. Curve 
a, Rayleigh's method using the mode for no ro- 
tation of a fixed-free blade, fundamental mode, 
/ 2 =2G.2 2 +1.85n 2 (c/s) 2 and second harmonic 
mode f /2=72.4H-6.45ft» (c/s) 2 ; curve b, Ray- 
leigh's method using the mode for no rotation 
of a free-free blade, fundamental mode,/ 2 =29.52 
+2.50n 2 (c/s) 2 , and second harmonic mode,/ 2 = 
83.2 2 +7.73n 2 (c/s) 2 . The closed points were 
values derived by Theodorsen from experi- 
ments with an affine model (one-tenth thick- 
ness) of blade type 4371. 



second harmonic frequency of pro- 
peller blades of type 32 and type 4371 
for the two extreme end conditions of 
rigid clamping (zero deflection, zero 
slope) at the hub and complete free- 
dom (zero slope, zero shear) at the 
hub. The results are listed in table 
1 and are shown as curves in figure 
16 (A) and 16 (B) for a blade of type 
32 and in figure 17 for a blade of 
type 4371. The open circles in the 
figures indicate the results of solving 
eq 23 and 24 for the natural mode. 
The solid points in figure 17 arethe 
results of experimental observations 
made by Theodorsen on a model 
blade of type 4371 [9], which results 
were communicated to the authors 
of the present paper by letter. 

The open circles for the "exact" 
solution are seen to check the curves 
calculated from Rayleigh's method 
within the computational error. The 
solid points obtained by Theodorsen 
lie inside the narrow band bounded 
by the curves for the fixed-free and 
the free-free end conditions. 
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Table 1. — Variation of natural frequency, f, with speed of rotation, n, for blades of 

two types 



Type of blade 


Method of deter- 
mining frequency 


End condition 


Mode 


Frequency, / 


32 


Rayleigb L__ _. 


Fixed-free 

Free-free 

Fixed-free 

Free- free. _ 


Fundamental 

do 

Second harmonic 2 ... 
do. 2 


c/s 

(30.92+1.95712)7 

(31.12+2.467*2)* 
(1 17.52+6. 8n2)* 

(1302+8.0n2)* 
(26.22+1.85712)* 
(29.52+2.50712)* 
(72.42+6.45712)* 
(83.22+7.737*2)* 

(302+2.05n2)* 
(84.12+6.07712;* 


32 


do.i 


32 


... do. 1 . 


82— 


do. 1 


4371 


do. 1 .. ... 


Fixed-free 

Free-free 

Fixed-free 

Free-free 


Fundamental 

do _._ 

Second harmonic 2 „_ 
do.2 


4371 


do. 1 


4371 

4371... 


do. 1 

do. 1 


4371 


From model blade 3 . 
do. 3 


Unknown 

.. do 


Fundamental 

Second harmonic 2 ... 


4371 









1 Using shape of natural mode for no rotation. 

2 Mode having one node near tip of blade. 

3 Test results by Theodorsen on model blade having one-tenth the thickness of prototype 4371. 

3. EFFECT OF RESTRAINT AT HUB; CORRECTION TO RP764 [1] 

Examination of table 1 and of figures 15 to 17 shows that a change in 
the condition of restraint of the blade at the hub may change the 
natural frequency of the flexural vibrations by as much as 13 percent. 
The natural frequency for the free-free end condition is 10 to 13 
percent higher than that for the fixed-free end condition for blades of 
both type 32 and 4371. 

An effect of restraint on frequency not considered in figures 15 to 
17 is that due to the finite length of the hub clamp. This will be 
equivalent to a shortening of the blade or to a displacement of the 
point of clamping toward the tip. Calculations for a nonrotating 
blade of type 32, whose point of clamping was moved as much as 15 
percent of the blade length toward the tip, gave an increase of 2 per- 
cent for the fundamental frequency as compared with the frequency 
for clamping at the center of the hub and an increase of 7 percent in the 
second harmonic frequency. The effect was too small to make it 
worth-while to extend the calculation to rotating blades. It appears, 
furthermore, that the direction of the effect is such as to increase the 
frequency for the clamped-free end condition so that the frequencies of 
blades with finite hub-clamps may be expected to lie inside the theo- 
retical bands given in figures 15 to 17. An experimental confirmation 
of this result is desirable, especially in view of the limitation of the 
theoretical work (see p. 644) to symmetrical modes with a restraint 
at the hub which may be replaced by a concentrated mass without 
any spring. 

The statement made above that the natural frequency of a pro- 
peller blade increases as the point of clamping is moved toward the 
tip, calls for a correction of results to the opposite effect cited on 
p. 205, 208, and 209 of [1]. An error was found in the calculation for the 
frequencies cited on these pages for a blade of type 32 clamped at a 
point removed from the center of the hub. Correcting this error 
showed that a displacement of the point of clamping equal to 15 per- 

100926—38 9 
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cent of the blade length raised the frequency for the fundamental from 
30.7 to 31.3 c/s and that of the second harmonic mode from 116.8 to 
125.0 c/s. A displacement of 5 percent in the point of clamping raised 
the frequency of the second harmonic mode from 116.8 to 121.8 c/s. 
After making these corrections it may be concluded that the resultant 
variation in frequency will lie inside the 13 percent range limited by the 
frequencies for the fixed-free and free-free end condition. This range 
is comparable to the variation of about 10 percent that was found in 
the maximum stress per unit tip deflection for different conditions of 
restraint. 

4. EFFECT OF INITIAL TWIST AND OF INITIAL CURVATURE 

The expression, eq 44, was derived for a straight propeller blade 
without initial twist or initial curvature. Actual propeller blades 
always show a large initial twist and a small amount of initial curva- 
ture. Keissner [17] concluded from his theory of twisted and curved 
propeller blades that the initial curvature had a small effect, while the 
initial twist may have a considerable effect on the natural frequencies 
of a rotating blade vibrating flexurally. He derived expressions cor- 
responding to eq 40, 41, and 42 for the potential energy and the 
kinetic energy of a twisted and curved blade vibrating flexurally while 
in rotation. These expressions show that a change in the angle of 
attack of the blade as a whole has no effect on p ] a twist varying along 
the blade has a certain effect. The smallness of the effect of initial 
twist on the fundamental mode and second harmonic mode of a non- 
rotating blade is indicated by the good agreement between the fre- 
quencies and stresses calculated for the straight blade and those 
measured for an actual twisted blade [1]. The smallness of the effect 
on the coefficient a occurring in eq 45 is shown by the close agreement 
between the calculated frequencies of the first two modes and those 
determined experimentally by Theodorsen on a model of a blade of 
type 4371. The principal axes of this blade twisted through approxi- 
mately 35° in passing from a section 8 inches from the center of the 
hub to a section 56 inches from the center of the hub. It may be 
concluded that the effect of twist on the fundamental and second 
harmonic frequency of a propeller blade of this design is small com- 
pared to the effect of end restraint at the hub. 

5. CRITICAL SPEEDS 

Curves of natural frequency versus speed of rotation, such as figures 
16 and 17, may be used to estimate the critical speeds of rotation of the 
propeller, that is, those speeds at which there may be resonance be- 
tween periodic impulses acting on the propeller with one of the natural 
frequencies of the propeller system. A condition of resonance is 
required to set up sufficiently high stresses in the propeller blades to 
cause their failure by fatigue. 

The frequency of the exciting forces acting on a propeller in flight 
will in general be proportional to the speed of rotation of the pro- 
peller. They will be represented by straight lines through the origin, 
such as are shown in figureg 18 and 19. The intersection of these 
straight lines with the bands denoting the natural frequencies gives 
the regions within which one may expect a critical speed. 
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The straight lines in figures 18 and 19 were drawn to illustrate the 
determination of critical speeds leading to flexural vibrations in the 
fundamental mode and the second harmonic mode for blades of types 
32 and 4371. 

Line a corresponds to exciting forces of the rotational frequency of 
the propeller, such as may be set up by a primary unbalance in the 
rotating system, or an asymmetry in the setting of one of the pro- 
peller blades which would cause a disturbance of propeller-rotation 
frequency in the airflow. This line does not intersect any of the blade- 
frequency curves since the blade frequency, /, is always higher than 
the frequency of rotation, n, due to the finite flexural rigidity of the 
blade. 

Line b corresponds to exciting forces having twice the rotational 
frequency of the propeller. Such excitation may exist in a two-blade 
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Figure 18. — Critical speeds for flexural vibrations with fundamental mode and 

second harmonic mode, blade type 32. 

Curve a, f=n; curve 6, /=2n; curve c, f='3n; curve d, f=3.5n; and curve e, f=7n. 

propeller in which each tip passes close to a part of the airplane struc- 
ture once per revolution. An exciting couple of twice the frequency 
of rotation will act on a two-blade propeller of an airplane flying in a 
curved path due to the gyroscopic couple [13] required to change the 
direction of the angular momentum of the propeller-engine combina- 
tion during the turn. 

Line c may enter in a three-blade propeller in which each tip passes 
once during a revolution close to a structural part of the airplane. 

Line d corresponds to the principal variation in the gas-pressure 
torque (3% times per revolution) for a single-row seven-cylinder 
engine with direct drive. 

Line e corresponds to the principal variation in the gas-pressure 
torque (7 times per revolution) for a double-row 14-cylinder engine 
with direct drive. 
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In geared engines one would have to multiply the frequency of the 
unbalanced torque components in the engine by the gear ratio to 
obtain the frequency of the torque pulsations acting on the propeller 
hub. 

The regions within which a critical speed may lie which leads to 
flexural vibrations in the fundamental mode or the second harmonic 
mode, are indicated by vertical shading. In the case of a blade of 
type 32, figure 18 indicates two dangerous regions for propellers rotat- 
ing between 1,000 and 2,000 rpm. Unbalanced torques from the 
individual cylinders in a double-row 14-cylinder engine may lead to a 
resonant tip vibration at a speed somewhere between 1,080 and 1,200 
rpm. Torque variations of the second order may excite resonant 
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Figure 19. — Critical speeds for flexural vibrations with fundamental mode and 

second harmonic mode, blade type 1+371. 

Curve a, f=n; curve b, f=2n; curve c, /=3n; curve d, /=3.5rc; and curve e, f=7n. 

vibrations in the fundamental mode somewhere between 1,340 and 
1,640 rpm. 

In the case of a blade of type 4371, figure 19 indicates two dangerous 
regions for propellers rotating between 1,000 and 2,000 rpm. | Eeso- 
nant vibrations in the fundamental mode may be excited by second 
order torque variations at a speed somewhere between 1,080 and 1,500 
rpm. Resonant tip vibrations may be excited by 3^-order torque 
variations due to the firing in the individual cylinders of a single-row 
seven-cylinder engine at a speed somewhere betweenll,800^and 2,300 
rpm. 

It should be emphasized that the above discussion is only intended 
as an illustrative example, and that it does not give a complete picture 
of all possible critical speeds that may be expected in a propeller of the 
type 32 and type 4371. Figures 18 and 19 do not contain the higher 
harmonics of flexural vibration which may be the cause of some tip 
failures in service [14]. Another omission is that mode in which the 
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propeller vibrates as a torsion pendulum with the crankshaft as elastic 
member. Its frequency in this mode may be calculated by well- 
established methods [15] ; it is nearly independent of the speed of rota- 
tion and will, therefore, be represented as a horizontal straight line in 
figures 18 and 19. Another mode that may be important in some cases 
is the torsional mode of the blade about its longitudinal axis. The 
frequency of this mode will be very high for all except very thin 
blades, such as solid-steel propeller blades. 

A further limitation which must be kept in mind in applying dia- 
grams of the type of figures 18 and 19 is that they neglect the coupling 
between different modes of vibration. The effect of coupling between 
the flexural vibrations of a propeller blade and the torsional vibration 
of the propeller about the crankshaft has been investigated for certain 
idealized cases by Carter [16]. 

In spite of all the limitations to which the diagrams in figures 18 and 
19 are subject, their introduction may serve a useful purpose in point- 
ing out the possibility of a large uncertainty in the critical speed cor- 
responding to a small uncertainty in the natural frequencies of a 
propeller. If all natural modes of vibrations were included in the 
figures, it is probable that only a small portion of the service range of 
speeds would lie outside the possibility of resonance. 

One will always have to contend with the presence of one or more 
critical speeds somewhere in the service range. It would be of some 
help if the location of these critical speeds could be determined much 
more accurately by experiment or otherwise than is possible with a 
chart such as figures 18 or 19, but it would be far better if the intensity 
of the periodic exciting forces acting on the propeller were reduced so 
that the stresses set up at a critical speed would lie below the endurance 
limit of the material. Adequate vibration insulation of the engine 
from the propeller would, for instance, reduce the excitation cor- 
responding to curves e and d in figures 18 and 19 to a low figure and 
may prevent fatigue failures due to tip vibration in the second har- 
monic mode. 

IV. COMPARISON OF NATURAL MODES FOR AFFINE 
PROPELLER BLADES 

The work involved in obtaining a solution for a natural flexural 
mode of a given propeller blade is too great for routine calculation. 
Complete solutions will probably be carried out for only a few charac- 
teristic types of propeller blades. This makes it important to find a 
model rule which would permit application of these particular solu- 
tions to propeller blades of other materials and of shapes that are 
affine 3 to the shape of one of the propeller blades that has been 
studied in detail. 

A model rule for predicting the natural frequencies of a given 
propeller blade from the frequency of a propeller blade of affine shape 
has already been stated and applied by Theodorsen [9]. Theodorsen 

s Two propeller blades, 1 and 2, are considered to bo affine in shape if the surface /t(ft, ft, z{) =0 of the first 
blade may be transformed to the surface fa (12, 1/2, ft)"«0 of the second by a linear transformation of the 
coordinates such that for corresponding points X2=aX\, y2=0yi, 2:2=7 t\, where a,0, 7 are constants. This 
may be expressed more concisely in mathematical terms by stating that the two surfaces may be written 
in the form/ (x, y, 2)=0and f (ax,0y r 7 z)=0, where a, 0, 7 are constants. 
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states that the frequency of a given mode of flexural vibration of 
similar propeller blades will vary as the ratio 

tjfc (48) 

where L is the length of the blade, D is a quantity representative of 
the thickness, say the diameter of the blade at the hub, and V B = (E/p) 1/2 
is the velocity of longitudinal sound vibrations in the blade material. 
For rotating propellers, Theodorsen derived the following model 
rule by proceeding from the frequency equation, eq 1. "The relative 
increase in frequency due to centrifugal force in similar 4 propellers 
is the same for the same value of the quantity 

SB- 

where V T is the tip velocity of the propeller." The use of the ap- 
proximate relation eq 1 as a proof for eq 49 is permissible in view of 
the excellent agreement between eq 1 and the more exact solution 
(see fig. 16 and 17). 

The model rules in eq 48 and 49 were derived for the flexural 
frequencies of propeller blades only. It is not safe, without further 
investigation, to conclude that they will give the conditions under 
which the flexural modes and hence the stress-distributions of two 
affine propellers are the same. An answer to this question for the 
idealized straight propeller blades, at least (see p. 641), may be 
obtained by determining under what conditions the deflection curves 
Xi{x x ) and X 2 (x 2 ) of two affine propellers satisfy the same differential 
equation, 8 or 10. If X 1 and X 2 satisfy the same equation at all 
points along the two propeller blades and if the conditions of restraint 
are the same, X x and X 2 must be proportional to each other according 
to eq 8 and the deflection curves and stress-distribution curves must 
be the same for the two blades. 

Let a, j8, 7 be the ratios of affinity for length, breadth, and thickness 
dimensions of the two affine propeller blades, that is, 

a=L 2 /L u P=b 2 lb h y=c 2 /c u (50) 

where L is the length of the blade and b and c are respectively the 
breadth and thickness of the blade at a section a given percentage 
of the blade length from the tip. The subscripts 1 and 2 refer to the 
blades 1 and 2 that are being compared. The condition that the 
deflection curves %=Xi/Li f and r) 2 =X 2 /L 2 of the two blades satisfy 
the same differential equation, eq 10, is expressed by 

^ 2 ^ 2 " - ^^+ (gl) J-^^^'] =^x^7x" - Ar^/^, + (^)^= ^> ']- (5 1) 

The ratios &, £ 2 > % h} > . . are, according to eq 9 and 50, related by 
fe-fu/^/i^-^i, & 2 =a 2 (g^)\, «,-£)k (52) 



* Theodorsen uses the word "similar" in [9] in place of the term "affine" used in this paper. 
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The substitution of eq 52 in eq 51 shows that eq 51 is satisfied if 

These model rules may be written in the more convenient form: 

The last equation is another statement of the model rule eq 48 that had 
already been given by Theodorsen for nonrotating propellers. Com- 
bining the second and third equations in eq 53, and noting that the 
tip speed, Vt=Lo), leads to Theodorsen's model rule, eq 49, for the 
relative rise in frequency due to rotation. 

It may be concluded that Theodorsen's model rules, eq 48 and 49, 
can be applied to calculate not only the frequency but also the shape 
of a natural flexural mode and hence the stress distribution of any 
idealized straight propeller blade from the corresponding natural mode 
of a "standard" propeller blade whose length, breadth, and thickness 
dimensions are affine to those of the propeller blade in question. 

Although Theodorsen's model rules were derived for idealized 
straight blades only, they are probably also applicable to actual 
twisted blades vibrating with low-frequency flexural modes, in view 
of the agreement found between the observed stress distribution in a 
nonrotating blade and that calculated for the corresponding idealized 
straight blade [1], and in view of the agreement found between the 
frequencies observed by Theodorsen on a model of blade type 4371 
and those calculated for the corresponding idealized straight blade 

An experimental investigation of the natural frequencies and the 
stress distributions in a rotating propeller blade vibrating in resonance 
with one of the flexural modes is highly desirable as a check on these 
conclusions as well as on the other theoretical results presented in this 
paper. 

V. SUMMARY 

The procedure for calculating the natural modes of flexural vibration 
of nonrotating propeller blades, which is described in RP764 has been 
extended to include the effect of rotation. The integral equations 
for the symmetrical modes of vibration of an idealized straight blade 
are given, in particular, for the cases of rigid clamping at the hub and 
zero clamping at the hub. These equations were solved numerically 
for aluminum-alloy propeller blades of two typical designs designated 
as type 32 and type 4371 rotating at speeds which cover the range of 
service speeds and vibrating with the fundamental flexural mode and 
with the second harmonic flexural mode (with a node near the tip). 

The effect of rotation in the operating range (1,000 to 2,000 rpm) 
on the maximum stress amplitude per unit tip deflection for the funda- 
mental mode was not more than 13 percent. It is sufficient for ap- 
proximate calculations to take the resultant stress in a rotating blade 
vibrating with a given tip amplitude as the sum of the static stress 
due to centrifugal force and steady air forces and a sinusoidally vary- 
ing stress whose amplitude is equal to that for the fundamental of the 
same blade vibrating without rotation with the given tip amplitude. 
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Application of the same procedure to the second harmonic mode 
leads to a still better approximation because of the decrease in the 
relative importance of centrifugal force at the higher flexural modes. 
There may be a shift towards the hub of as much as 20 percent of the 
blade length in the location of the point of maximum stress amplitude 
in those cases in which the stress for no rotation decreases slowly on 
the hub side of the point of maximum stress amplitude. 

The natural frequencies of a rotating propeller blade vibrating 
flexurally can be calculated by solving the integral equations. This 
is usually too cumbersome for practice and it is unnecessary since an 
excellent and convenient approximation to the frequencies obtained 
in this manner follows from Rayleigh's method by substituting the 
deflection curve of the propeller blade for no rotation in the expressions 
for potential and kinetic energy. Rayleigh's method was used to 
calculate the natural frequencies of the fundamental mode and the 
second harmonic mode for both types of propeller blades and for the 
two end conditions of rigid clamping and no clamping at the center 
of the hub. The frequency of the two flexural modes was found to 
change 9 to 13 percent in passing from no clamping to rigid clamping. 
Observed values of natural frequency for the fundamental mode and 
the second harmonic mode of a model blade of type 4371 which were 
obtained by Theodorsen at the Langley Memorial Laboratory were 
found to lie inside the narrow band bounded by the curves for rigid 
clamping and no clamping. 

A displacement of the point of clamping from the center of the hub 
toward the tip such as may be expected due to the finite length of hub 
clamp will slightly increase the natural frequency of flexural vibra- 
tions. This calls for a correction to RP764 [1], where the opposite 
conclusion had been reached due to computational error. 

The natural frequencies determined by Rayleigh's method may be 
used to calculate the critical speeds for flexural vibration. Unfortu- 
nately, it appears that a given error in the natural frequency of the 
blade may introduce several times that error in the critical speed, so 
that the calculated regions containing a critical speed may cover a 
considerable range of speeds. 

Both the natural frequencies for flexural vibration and the stress 
distribution of a given idealized straight propeller blade may be ob- 
tained from the observed frequencies and stress distributions of a 
model blade that is affine to the blade in question by making use of 
the two model relations which were first applied by Theodorsen to 
calculate the frequencies of a propeller blade of type 4371 from a 
model blade having the same length and width dimensions but having 
one-tenth the thickness of the prototype. 

The above results are based upon theoretical reasoning and hence 
upon certain idealizing assumptions that can only be realized approx- 
imately in an actual propeller blade in flight. An experimental in- 
vestigation of the natural frequencies and the stress distributions in a 
rotating propeller blade vibrating in resonance with one of the flexural 
modes is highly desirable as a check on these predictions of the theory 
and to establish the modes of vibration which are actually important 
in setting up severe alternating stresses in a propeller. More experi- 
ments are needed to determine the effects of twist on the natural 
flexural modes; existing tests have shown this effect to be small only 
for two nonrotating blades and for the fundamental and second 
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harmonic frequencies of one rotating model blade. Experiment only 
can give the actual amplitudes of the resonant vibrations; the theory 
may be of assistance in deducing the amplitude of stresses throughout 
the blade from the vibration amplitude measured at a convenient 
point. The carbon resistance strain pick-ups recently developed by 
A. V. deForest at the Massachusetts Institute of Technology may 
make an adequate experimental study of propeller vibration possible 
in the near future. 
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